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I. INTRODUCTION 
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We present a detailed calculation of our previous short paper [M. Shiraishi, D. Nitta, S. Yokoyama, 
K. Ichiki, and K. Takahashi, Phys. Rev. D 82, 121302 (2010).] in which we have investigated a 
constraint on the magnetic field strength through comic microwave background temperature bispec- 
trum of vector modes induced from primordial magnetic fields. By taking into account full angular 
dependence of the bispectrum with spin spherical harmonics and Wigner symbols, we explicitly show 
that the cosmic microwave background bispectrum induced from the statistical-isotropic primordial 
vector fluctuations can be also described as an angle-averaged form in the rotationally invariant 
' way. We also study the cases with different spectral indices of the power spectrum of the primordial 

\ magnetic fields. 

^ ■ PACS numbers: 98.80.Cq, 98. 62. En, 98.70.Vc 

o 

\ Recent observational consequences have shown the existence of 0(1O~ 6 )G magnetic fields in galaxies and clusters 
of galaxies at redshift z ~ 0.7 — 2.0 [H-Q- One of the scenarios to realize this is an amplification of the magnetic fields 
Oh by the galactic dynamo mechanism (e.g. Q), which requires 0(1O~ 2O )G seed fields prior to the galaxy formation. 
A variety of studies have suggested the possibility of generating the seed fields at the inflationary epoch p| Q , the 
cosmic phase transitions 0, H| , and cosmological recombination and also there have been many studies about 

the constraint on the strength of primordial magnetic fields (PMFs) through the impact on the cosmic microwave 
background (CMB) anisotropics, in particular, the CMB power spectrum sourced from the PMFs [T34I3- 

Recently, in Refs. [18-22], the authors investigated the contribution to the bispectrum of the CMB temperature 
CNj ■ fluctuations from the scalar mode PMFs and roughly estimated the limit on the amplitude of the PMFs. Because 
the temperature fluctuations induced by the PMFs have the highly non-Gaussian statistics, the bispectrum of such 
fluctuations should have nonzero value. As is well known, PMFs excite not only the scalar fluctuation but also the 
vector and tensor fluctuations. In particular, it has been known that the vector contribution may dominate over 
the scalar one on small scales by the Doppler effect in the CMB power spectrum (e.g. [12, [HI). Hence, the future 
CMB experiments, for example, Planck satellite [23|, are expected to give a tighter constraint on the amplitude of 
the PMFs from the vector contribution induced from the magnetic fields. With this motivation, in Ref. [24|, we have 
presented a CMB angle-averaged bispectrum of the vector perturbations induced from the PMFs and also a forecast 
of upper limit for the strength of PMFs smoothed on IMpc scale as -BiMpc < 10nG. However, there we could not 
• • , show the details of calculation for the limit of pages. Hence, in this paper, we focus on the derivation of the CMB 
bispectrum of vector modes induced from PMFs without neglecting the full angular dependence on the wave number 
vectors. In this paper, we also show that the CMB vector bispectrum induced from the statistical-isotropic PMFs 
can be described as the angle-averaged form like the scalar mode (25j in the rotationally invariant way 1 . 

This paper is organized as follows. In the next section, we formulate the CMB vector bispectrum induced from 
PMFs by following the procedure of Ref. [Hj]. In Sec. IIII1 we analytically expand the CMB bispectrum with help 
from some numerical evaluations. In Sec. IIV1 we show our result of the CMB bispectrum from the PMFs and estimate 
the limit of the amplitude of the magnetic fields. In addition, we also discuss the shape of the bispectrum. The final 
section is devoted to summary and discussion of this paper. 

Through this paper, we assume the universe is spatially flat and use the definition of Fourier transformation: 

d 3 k 



/<7 J k 
(2^ /(k)eJk X ' (1) 



'Electronic address: marc@a. phys.nagoya-u.ac.jp 

1 In Ref. |2(J , the authors presented the analytical formulas of the CMB vector bispectrum sourced from statistically isotropic PMFs in 
a different approach than ours and claimed that the bispectrum violates the rotational invariance. Recently, however, they also could 
reduce the final formulas to the rotational-invariant form, which will be shown in an updated version of their paper |27H . 
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II. FORMULATION OF THE VECTOR BISPECTRUM INDUCED FROM PMFS 



Let us consider the stochastic PMFs B b (x, r) on the homogeneous background Universe which is characterized by 
the Friedmann-Robertson- Walker metric, 



ds 2 = a(r) 2 [-dr 2 + 5 bc dx b dx c ] 



(2) 



where r is a conformal time and <z(t) is a scale factor. The expansion of the Universe makes the amplitude of the 
magnetic fields decay as 1/a 2 and hence we can draw off the time dependence as B b (x, r) = B b (x)/a 2 . The Fourier 
components of the spatial parts of the PMFs' energy momentum tensor (EMT) are described as 



Tt(k,r)Ep 7 (r)[^A fl (k)+4 c (k)] , 
1 f d 3 k' 



A B (k) = 



87rp 7i o 



n fc Bc (k) = - 



i 



(2tt) 3 
d 3 k' 



B»(k')B b (k-k') 



(3) 



47rp 7i0 J (27r) 



B b (k')B c (k~k') , 



where we have introduced the photon energy density p 7 in order to include the time dependence of a and p 7 ,o 
denotes the present energy density of photons. In the following discussion, for simplicity of calculation, we ignore the 
trivial time-dependence. Hence, the index is lowered by S bc and the summation is implied for repeated indices. 
Assuming that B a (x) is a Gaussian field, the statistically isotropic power spectrum of the PMFs Pb(/c) is defined 

by 



(B a (k)B b (p)) = {2nf^P- 



P Qb (k)5(k + p) , 



(4) 



with a projection tensor 



P ab {k) = V ei°Hk)e [ - a) (k) = S ab - k a k h 



cr=±l 



(5) 



which comes from the divergence free nature of the PMFs. Here k denotes a unit vector and ei is a normalized 
divergenceless polarization vector which satisfies the orthogonal condition; ftj 1 ' = 0. The details of the relations 
and conventions of the polarization vector are described in the Appendix in our previous paper (28| . Although the 
form of the power spectrum Ps(k) is strongly dependent on the production mechanism, we assume a simple power 
law shape given by 



P B {k) = A B k T 



(6) 



where Ab and ns denote the amplitude and the spectral index of the power spectrum of magnetic fields, respectively. 
In order to parametrize the strength of PMFs, we smooth the magnetic fields with a conventional Gaussian filter on 
a comoving scale r: 



B. 



k 2 dk 
~2^ ( 



P B {k) , 



(7) 



then, Ab is calculated as 



An = 



(27T) 



n B +5 



T(H£±^)kl'! 



(8) 



where r(x) is the Gamma function and k r = 2ir/r. 

We focus on the vector contribution induced from the PMFs, which comes from the anisotropic stress of the EMT, 
i.e., Tisab- Using the polarization vector, the vector anisotropic stress fluctuation is given by 



H^ 1) (k) = fc a 6i Tl) (k)H Bab (k). 



(9) 



In the magnetic case, this acts as a source of the CMB fluctuations of vector modes. 
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A. Bispectrum of the vector anisotropic stress fluctuations 

As we have mentioned above, the PMF B b is assumed to have Gaussian statistics. Hence one can easily find that 
the statistics of the vector anisotropic stress fluctuation given by Eq. (0) are highly non-Gaussian and the bispectrum 
(3-point function) of that has a finite value. 

Using Eq. (0} and the Wick's theorem, the bispectrum of n^ 1 ^ (k) is calculated as 



(n BQ f,(k 1 )n Bcd (k 2 )n Be/ (k 3 )) = (-^p 7 , ) 



n Sa6 (k 1 )n Bc(i (k 2 )n Be/ (k3)) fc lae ^ Al, (k 1 )fc 2c e^ A2, (k 2 )fc 3e e> A3) (k 3 ) , (10) 



n 



n 



(2tt) 3 



Ln=l ' 

x (B a (ki)B 6 (ki - k£)B c (ki)B d (k a - k^B e (k' 3 )B f {k 3 - k' 3 )) 

!.'2 Jul r>_nJ \ I j2i 



(-4vrp 7! o) 



-3 



n 

.n=l 



k'^dk' n P B {k' n ) / d 2 k' n 







x<5(k! - ki + k 3 )<5(k 2 - k' 2 + ki)5(k 3 - k 3 + k' 2 ) 

x l[P ad (k\)P be {k' 3 )P cf (k\) + {a o b or c o rf or e o /}] 



(11) 



where A„ denotes the helicity of the vector mode as A„ = ±1 and Ud is the Alfven-wave damping length scale [23. [301] 
as A;^ 1 ~ O(0.1)Mpc and the curly bracket denotes the symmetric 7 terms under the permutations of indices: a o b, 
c H ii, or e (4 /. Note that we express in a more symmetric form than that of Ref. fl8| to perform the angular 
integrals which will be described in Sec. IIIII To avoid the divergence of (n,B a ft(ki)nB C d(k 2 )nB e /(k 3 )) in the IR limit, 
the value range of the spectral index is limited as iib > —3. 



B. CMB all-sky bispectrum 

The CMB temperature and polarization fluctuations are expanded into (spin-weighted) spherical harmonics (l2l 

-,(Z) 



1311 . |32[ . Then the angle-averaged bispectrum formed by their coefficient a x tm can be defined as [25j, |2£ 

/ 3 \ 



B 



X 1 X 2 X 3 ,£ 



E 



ii 12 h 

mi 1112 ^3 



\n=l 



(Z») 
l X n ,i„m„ 



(12) 



where the matrix is the Wigner-3j symbol, Z = S, V or T is corresponding to the scalar, vector or tensor-mode 
perturbation respectively, and X = I, E or B means intensity, .E-mode or i3-mode polarization, respectively. 

Let us consider B^Yt^i^ induced from 11^. In the same manner as in Refs. [HI, [Hj], ci\} m sourced from PMF is 
given by 



i {V) 



47r(— i) 



A=±l 



(2^) 3 



n^imW = / d 2 k4 ± „ 1) (k) Tl F; m (k) 



(13) 
(14) 



Here T^P denotes the radiation transfer function of the temperature fluctuation from magnetic vector mode as 
calculated in Appendix \K\ and ^iY« m (k) is the spin-1 spherical harmonic function. By making use of these equations, 
we obtain the CMB temperature bispectrum induced from the vector anisotropic stress n^ 1 ^ which is given by 



R (VVV) _ 



T^(k 1 ,k 2 ,k 3 )^(2n)- 3 



ti h e 3 

mi m,2 TO3 



nn 



Ai,Aj,A 3 =±l 



(15) 
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These equations are corresponding to Eqs. (2.9) and (2.7) of Ref. [Hj]. 

In the next section, we will derive an explicit form of B^^^J 2i by calculating the complicated angular dependencies 
on the wave number vectors, which are implied by Eqs. (|10p . (|14[) and (|15[) . with the spin- weighted spherical harmonics 
and the Wigner symbols. In the calculation, we will see that the dependence on the azimuthal quantum numbers 
(mi, m 2 and m 3 ) in the bispectrum of II^ j 1 ] m is confined only in the same form as the Wigner-3j symbol in Eq. (|12l) . 
which implies the rotational invariance of the CMB bispectrum from the vector anisotropic stress of PMFs [24[ . 



III. ANALYTIC CALCULATION OF THE CMB TEMPERATURE BISPECTRUM 

In this section, we derive the explicit form of Eq. (fT5j) by calculating the full-angular dependence which has never 
been considered in the previous studies fl9T - l22l f26j j . The following procedures are based on the calculation rules 
discussed in Ref. [28} . Note that we use some colors in the following equations for readers to follow the complex 
equations more easily. 



A. Exact expression of J~^\ 3 3 

Let us consider an exact expression of ^ | 3 , by expanding all the angular dependencies with spin- weighted 
spherical harmonics and rewriting the angular integrals with summations of angular and/or azimuthal quantum 
numbers. Substituting the expression of the bispectrum of Il^^k) (Eq. ([TTJ)) and Eq. flU]) into Eq. §T5§, we can 
obtain 



E 



«1 «2 «3 
TOi m 2 TO3 



xj(k! - K + K)S(k 2 - k' 2 + k' ± )6(k 3 k 3 + k / 2 )k\J- Xl \k\)k 2 J d ' X2 \k\)k 3 Jf X3 \k 3 ) 



P ad (K)P be (k' 3 )P cf (k' 2 ) + {a <-> b or c <+ d or e <+ /} (-87^,0)" 



(16) 



At first, we focus on the first term of permutations. 

In the first step, in order to perform all angular integrals, we expand each function of the wave number vector with 
the spin-weighted spherical harmonics. By this concept, three delta functions are rewritten as [H, l33j 



<J(ki - ki + k 3 ) = 8 / A 2 dA V (-1 

JO T..T.T.. 



IIXlJl, (hA)j L2 (k[A) jL3 (k' 3 A) 



S(k 3 



L1L2L3 
M1M2M3 



<y(k a -ki + ki) = 8 



B 2 dB £ (-l) L ' 1+ ^ +L ^ . 



l:LlL^L'Sk2B) 3L ,(k' 2 B)j L ,{k[B) 

1 2 3 

x Y£, M , (k a )y £ , M , (k' 2 )Y£, M , (k\)(-l) M -- 



L-^ -^2 ^3 
M{ -M 2 M3 



/*°° L" +3L /f + h / ' 

k' 3 + k' 2 ) = 8 C 2 dC ]T (-1) 1 3 ' 3 ^L»JLy(fc 3 C')j L »(^C')j L »(fc^) 



x^ M „(k3)^M^(k 3 )^ M ,,(k 2 )(^i)^f ^, 



(17) 



where 



TS1S2S3 



(2^ + 1) (2^ 2 + l) (24 + 1) / ^ 4 4 
47r \ si S2 S3 



(18) 
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The other functions in Eq. (|16l) . which depend on the angles of the wave number k„, can be also expanded in terms 
of the spin-weighted spherical harmonics as 

CT=±1 

= E E (y) (-A 2 )r lma (k;)_ A2 r lmd (k 2 )_ ff r 1 ;„jki) CT y 1 ;„ d (ki) , 

/ 4 \ 2 ( 19 ) 

fc 2c e^ A3) (k 3 )p c/ (4) = 2 E hr (-A3)yi mc (k2)-A3n m/ (k 3 )^r 1 * m jk^) CT ,Y 1 ; n/ (k^ 



<t"=±1 m c ,m,f 



1) > 



4tt x 2 



fc 3e 4 Al) (k 1 )P 6 e(k 3 ) = ^ ^ (t) (- A i) y ^( k 3)-^ y i^(kO-'nk(kl) CT ^i* mt (kl) , 

cr'— ±1 m(,,m e ^ ' 

where we have used Eq. ([5]) and the notations of a unit vector n and a divergenceless unit vector e£ as [111 

"a = ^a™Yi m (n) , 

771 

^(n) = 4^>*(n) = T £ a™ ±1 Y lm (n) , ( 20 ) 

771 

/ 47T 

„ m m / i \ m c 

"a "a = °» l ' 

In the second step, let us consider performing all angular integrals and replacing them with the Wigncr-3j symbols. 
Three angular integrals with respect to k ^ , k' 2 and k^ are given as 

o-a-s 7-0-O--S ( L' 3 1 L \ ( L 2 1 L 



f d 2 k\^Y* m Y L2M2 ,Y* md Y*, M ,=J2 E (- i r +maI ^iL Sl L 2 : 

LM S=±l 



Mi m d M / \ M 2 -m a M I ' 



Jd 2 k\- a ,,Y* m Y L , 2M , a ,,Y* mf YZ,, M ,,- ^ ^ ( 1) CT +'" =/°„^, 5 m/ M ,^J _ mc M ,^J ,(21) 



L'M' S'=±l 

ii i r " 



I d*k\^Y? m Y L ,, M ^Y? mb Yi 3M3 = e E (-i) CT ' +m ^U- s "^rr s " 1 Lz L " ^ ( L " 

•J T II A III Oil Li 



ii" I M 3 m b M" / I Ml! -m e M' 

L"M" S"=±l 



where we have used a property of spin- weighted spherical harmonics given by [28l [3j 

' h h h 
mi m 2 to 3 



n^ Fz " m " - E •s^ismaVfe h "( m \ mo rL I ' ^ 22) 

71 — 1 /37773S3 

YL = (-i) s+m _ s y ; _ rn . (23) 



s ± lm 

We can also perform the angular integrals with respect to ki , k 2 and k 3 as 



/ d 2 kx-\ 1 Yi mb Yx ma _x- l Yl imi Y^ Mi —EE ^L^Lk^- 

" T AT O Li 



11L, 

L k M k S k =±l 



rP-V V V v* V* - r 0X 2 ~S p T 0X 2 -S p I L,\ 1-2 Li p \ I 1 1 L p \ 

a k 2 _ X2 r lmd Y lmc _ X2 Y e2m2 Y L , iM , - 2^ 2^ h'^L.hiL, y m[ m 2 M p J \m r m d M n (24) 

rP-H V V v* v* r 0A3-S 9 r 0A 3 -S, / ^ 3 L q \ f 1 1 Lg 

a K 3 ^ X3 Y lmf r lmc _ X3 Y l3m3 Y L ,, M ,, - 2^ 2^ L viHL q hiL q [ M n m Mn ) \ m F mf M n 

L a M a S a =±l 



Li 


h 




Mi 


mi 


M k 


L'i 


(2 


Lp 


M{ 


m 2 


Mp 


L'l 




L q 


M'{ 


m 3 


M q 



At this point, all the angular integrals in Eq. ([16]) have been reduced into the Wigner-3j symbols. 

Then, in the third step, we consider summing up the Wigner-3j symbols in terms of the azimuthal quantum numbers 
and replacing them with the Wigncr-6j and 9j symbols, which denote Clcbsch-Gordan coefficients between two other 
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eigenstates coupled to three and four individual momenta [281. |34M36| . Using these properties, we can express the 
summation of five Wigner-3j symbols with a Wigner-9j symbol: 

Li L 2 L 3 \ f 1 1 L k \ ( L 3 1 L" \ ( L 2 1 L \ f L 1 £ 1 L, 

k 



MiM 2 m 3 



\ Ah -M 2 M 3 M m„ m 6 M fc J I M 3 m b M" ) I M 2 -m a M J I Afi mi M, 



Al p m c md 




A4+,n r f i'l ^2 ^3 V 1 1 L p\[ L 3 1 L \( L 2 1 \ ( L[ l 2 L p 

[ > \ M[ —M2 M3 y ^ m c m d M p y ^ M3 m d M y y M^ -m c M' J \ M[ m 2 M p 



M'+e 2 +L'+L' + l 



M['M^'M 3 
M q m e mf 




(25) 



= (-1) 

V (_^+m e ( L'l L 2 ' L'i \( 1 1 L q \( L'i I L'\( L' 2 > 1 L" \ / L" l 3 L q 
2.^ 1 LJ \ M'( -Mil M'l ) \ m e m f M q I I M'l m f M' J I M'l -m e M" ) \ Ml m 3 M q 



(-1) 



M"+e 3 +L"+L" + l 



L' L" l 3 
M' —M" m 3 




Furthermore, we can also sum up the renewed Wigncr-3j symbols arising in the above equations over M , M' and M" 
with the Wigner-6j symbol as [37| 



L" L A / L V i 2 \ f L' L" l 3 
M" -M mi ) \ M —M' m 2 ) \ M' -M" m 3 



_^\M+M'+M' 

' h i 2 l 3 \\ h h I 



' \ mi m 2 m 3 / I L L L J 

With this prescription, one can find that the three azimuthal numbers are confined only in the Wigncr-3j symbol as 

£1 £ 2 £ 3 \ , . . . (±1) 

. This 3? symbol arises from the bispcctrum of II V, i and exactly ensures the rotational invariancc 
m 1 m 2 m 3 J J J 1 Bv,tm •> 

of the CMB bispectrum as pointed out above. 

So far, we have considered only the first term of permutations in Eq. (|16p . Hence, finally, we have to consider the 

contribution of the other 7 permutations. For example, in the calculation of the {a 6} part, m a and mt> of Eq. (|2ip 

replace each other and the summation over m a and mj, in Eq. (|25[) changes as 

Li L 2 L 3 \ ( 1 1 L k \ f L 3 1 L" \ ( L 2 1 L \ ( L x £ x L 



JUi M 2 M 3 
M k m b m a 



ST (_i\M 2 +m b I ^1 ^2 ^3 W 1 I I I / 1 ^ W 1 1 fc 

I Mi -M 2 M 3 M m m 6 M fe J I M 3 m a M" J { M 2 -m b M J I Mi mi M fc 



M) M+£l+i3+L+1+ifc ( ^ _w ™. ) < ^ L 2 ii > , (27) 




hence, the extra factor (— l) Lk arises. In the same manner, we can find the extra factor (— l) Lp or (— l) Lq in the 
{c <H> d} or the {e O /} part, respectively. 

Using the above expansions and the orthogonality of the Wigner-3j symbols given by 



/ «1 «2 «3 



01111121113 



mi m 2 m 3 

and performing the summations over Lk , L p and L g such as 



= 1 , (28) 





ECl^- Al± ^<! L 3 L 2 Li )=--^=ll% X { ^ L 2 Li\, (29) 
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we can obtain an exact form of J 7 ^^ 3 given by 

x E E 

where 



LL'L" S,S',S"=±l 



1 h h \ fS"SX lf ul V U.\fSS'X 2 fu' V u \fS'S"Xaf U i U l 



L L" L 



fu'Lil (^3 1 M.) k\)f 1,1,1 £ {k'i, fc 2 , L i L ng 3 (fc 2 , k 3 , ks){30) 



2(8tt) 



3/2 



/>oo 

V / A 2 dAj i 3(r 3 A)j i2 (r 2 A)i il (r 1 A) 

T T JO 



L1L2L 



(31) 



1 1 2 



This expression is one of our results in this paper. The above analytic expression of Fg^i^ 3 seems to be quite useful to 
calculate the CMB bispectrum of vector modes induced from PMFs with the full angular dependence. However, it is 
still too hard to calculate numerically, because the full expression of the bispectrum has six integrals and summations 
over the helicities as 



B 



(VVV) 
III Mils 



8(2tt) 



1/2 



E 

LL'L" 



tl h h 

V L" L 



E 



L1L2L3 



L" L l x I [ L L 1 i 2 

1 LilIl Jl>%Il'Jl>{l\l>; { Lz L2 L } \ { L } L 2 L l } \ L 3 L 2 L l 



1 1 2 



1 1 2 



3 «QO 7,2 JL poo pOO pOO 

/ TTriT^Hh) / A'dAj^ihA) / B 2 dBj L{ (k 2 B) / C 2 dCj L ,,(k 3 C) 
i=1 Jo ( Z7T ) J Jo Jo Jo 

x / ° k[ 2 dk[P B (k[)j L2 (k[A)jLdk[B) ' k?dk 2 P B (k 2 )j L ,MB)j L ,,(k' 2 C) 



x / kfdk' 3 P B (k' 3 )j L ,,(k' 3 C)j L3 (k' 3 A) 

E/_1 \L 2 +L 2 +L'1+L S +L' 3 +L'3 tOS~S tOS-S jOS'-S' tQS'-S' tOS"-S" tOS"-S" 
^ i / iUi i-2li 1 L'.'1L< L' 1L' 1 L 3 1L" 1 L"IL" 



S,S'S"=±1 



EiflArA, rOA 2 -A 2 7-OA3-A3 
i Li«i2 J, L' 1 £ 2 2 i L' 1 '^ 3 2 



(32) 



Ai,A 2 ,A 3 =±l 



In the following subsection, we introduce an approximation, the so-called, thin last scattering surface (LSS) ap- 
proximation to reduce the integrals and perform the summations over the helicities based on the selection rules for 
Wigncr-3j symbols. 



B. Thin LSS approximation 

Let us consider the parts of the integrals with respect to A, B, C, k' ,p' and q' in the full expression of the bispec- 
trum (Eq. (|32|l ). In the computation of the CMB bispectrum, the integral in terms of fc (, p and q) appears in the 
form as J k 2 dkTj^(k)jL 1 (kA). We find that this integral is sharply-peaked at A ~ tq — r*, where To is the present 
conformal time and r* is the conformal time of the recombination epoch. According to Ref. fl2l [l3[ , the vorticity of 
subhorizon scale sourced by magnetic fields around the recombination epoch mostly contributes to generate the CMB 
vector perturbation. On the other hand, since the vector mode in the metric decays after neutrino decoupling, the 
integrated Sachs- Wolfe effect after recombination is not observable. Such a behavior of the transfer function would 
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FIG. 1: (color online). The ratio of the left-hand side (exact solution) to the right-hand side (approximate solution) in Eq. 
(|33|l . The lines correspond to the case for L\ — l\ + 2 (red solid line), for L\ — l\ (green dashed line), and for L\ = l\ — 2 
(blue dotted line). 



be understood based on the calculation in Appendix lAl and we expect TpQ(k) oc ji x {k{ra — t*)), and the fc-integral 
behaves like S(A — (to — 7*)). By the numerical computation, we found that 



Ad A 



kfdhTj^ih^ihA) ~ (r - r») 2 f kfdhT^ih^eAhiro ~ 7-0) , 



(33) 



is a good approximation for L\ — £% ± 2, £i as described in Fig. [TJ Note that only the cases L\ = £i ± 2, i?i should be 
considered due to the selection rules for Wigner-3j symbols as we shall see later. From this figure, we can find that the 
approximation (the right-handed term of Eq. (f33|)) has less than 20% uncertainty for l\ ~ L\ > 100, and therefore 
this approximation leads to only less than 10% uncertainty in the bound on the strength of PMFs if we place the 
constraint from the bispectrum data at £1,^2^3 ^ 100 2 . Using this approximation, namely A = B = C — > To — t„ 
and J dA = J dB = J dC — > t»/5, the integrals with respect to A, B, C, k' ,p' and q 1 are estimated as 

3 coo .2, 



k n dk n rr (y} /u . 

(2tt)3 'i^ Kn > 



A dAji J1 (kiA) / B 2 dBj Lli (k 2 B) / C 2 dCj L ,,(k 3 C) 



[° kfdk'.PBik'^jLM^jL^KB) f° k' 2 2 dk' 2 P B (k' 2 )3 L ,(k' 2 B) jL ,,(k' 2 C) 
Jo Jo 



k' 2 dk' 3 P B (k' 3 )j L; ,(k' 3 C)j L3 (k' 3 A) 



~i 1 (k n )je n (k n {T a - r*)) 



(to-t*) (^yj /C ia ^, '(to-t*)/^,, ; (r - T,)lC L h L3 ; (r -r,) 



(34) 



Here the function JC,,, is defined as 



7T ?/ 



dfcfc 

JV-l 



l-N 



ji{ky)ji'{ky) 



T(N)T( l+l '+ 2 2 - N ) 



2y 2 N T( l ~ l '+ 1 + N )r( - W+i+w )r( W+2+w ^ 



(for y, TV, Z + /' + 2 - TV > 0) 



(35) 



2 Of course, if we calculate the bispectrum at smaller multipoles, we may perform the full integration without this approximation 
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which behaves asymptotically as JC^ v (y) oc I N for I ~ I' ^> 1. Here we have evaluated the k' integrals by setting 
kr> — > oo. This is also a good approximation because the integrands are suppressed enough for k',p',q' < kjj ~ 
O(10)Mpc- x . 



C. Selection rules of the Wigner-3j symbol 



Next we consider performing the summations with respect to the helicities of vector modes. By considering the 
selection rules of the Wigner-3j symbol, the summations over S, S' and S" (red part in Eq. (|32[1 ) are performed as 



E(_1 \L2+L 2 +L%+L 3 +L' 3 +L 3 ' tOS-S tOS-S tOS'-S' tOS'-S' tOS"-S" tOS"-S" 
y 1 ) ' L' 3 1L 1 L 2 IL L'^IL' L' 2 1L> 1 L 3 1L" 1 U 2 '\L" 

S,S',S"=±1 

= f^r^ll^^&^li" for L 3 + £2,^' + ^3 + ^' = even 
[ otherwise 

By the same token, the summations over Ai, A2 and A3 (green part in Eq. (|32p) are given by 



(36) 



\p / oa 1 -a 1/ oa 2 -a 2/ oa 3 -a 3 = i ^lle^v^ 1 !} e 3 \ for L i + ^ L 'l + **> L l + ^ = even 
, ^ ^ Lll!l2 L ' lt22 L "* 3% |0 otherwise ' 1 

Ai,A2,A3 — ±1 V. 

Then, using the function JCff, and the above equations, the CMB bispectrum of Eq. (ITS]) can be written as 

"°° 1,2, 



<wv) „ ( Z2{2k) 1,2 \ 
D III,Mil 3 — 



3/0 7 ,o / 



II 4 - [ ^T^ik^uMro - r.)) 



v rOl-l rOl-1 rOl-1 

X 1 L 1 t 1 2 1 L' 1 £ 2 2 1 L';i 3 2 

L\L'-^L'-^ LL' L 

3 



4 4 ^3 



LqLq l$ 

1 \E 3 -i ii±£i^i±£t r o j-o 7-000 7-01-1 7-01-1 7-01-1 7-01-1 7-01-1 7-01-1 

X l " 1 L 1 L2L3 1 L' 1 L' 2 L'/L'{L^L^ 1 L'lL 1 L2lL 1 L''lL' 1 L' 2 lL' 1 L3lL'' 1 L^lL'- 




(38) 

Here from the selection rules of the Wigner symbols [28j , we can further limit the summation range of the multipolcs 

as 

\h-h\<^<h+^ , 

Li = \£i±2\, l x , Li = |i 2 ±2|, 4, L? = |4»±2|, 4» , 

|£-4| < L' < i + , Max[|£-4|,|^'-4|] <i"<Min[L + 4,i' + 4] , 

(£ 2 ,L 3 ) = (|L-1|,|L±1|), (L,L), (L+1,|L±1|), 

(^,^) = (|L , -1|,|L'±1|), (I/,Z/), (i' + l,|£'±l|) , 1 j 

(L' 2 ',L 3 ) = (|L"-1|,|L"±1|), (i",i"), (L" + 1,\L"±1\) , 

Li + L 2 + £3 = even , L\ + L' 2 + L' 3 = even , L'[ + L 2 + L 3 = even , 

\L 1 -L 2 \<L 3 <L 1 + L 2 , \L[ ~ L' 2 \ < L' 3 < L[ + L' 2 , \L" — L 2 \ < L' 3 ' < L" + L 2 . 

and from the above restrictions the multipoles in the bispectrum, £\,£ 2 and £3, are also limited as 

£1 + £2 + £3 = even . (40) 

Therefore, these selection rules significantly reduce the number of calculation. In these ranges, while L' and L" 
are limited by L, only L has no upper bound. However, we can show that the summation of L is suppressed at 
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4 ~ 4 ~ 4 <C L as follows. When the summations with respect to L, L' and L" are evaluated at large L, V and L", 
namely 4,4, £3 -C L ~ L' ~ L", L 2 ~ L' 3 ~ L, L' 2 ~ L 3 ' - L' and L 2 ' ~ L 3 ~ L", we get 



EJ 4 4 4 \ r-("s +1 )^ ^. \r-("B+ 1 )/, 



LL'L" K J L 2 L' 2 L 2 



\ <, — (n B +l) / \ 




, r_ -nZ) 3 = i £ " i+ - L j +L ^' ,000 7-000 7-000 7-01-1 7-01-1 7-01-1 7-01-1 7-01-1 7-01-1 

' v 1 ) 1 L 1 L^L 3 1 L' 1 L' 2 L' 3 1 L'{LqL'.i 1 L' 3 lL 1 L 2 lL 1 L^lL' 1 L' 2 lL' 1 L 3 lL" 1 L^lL" 

L L' £2 I f L' L" 4 

^3 ^2 L[ > < L3 I/2 L" 

1 1 2 J [ 1 1 2 

oc ^ (LL'L") nB+4/3 ■ (41) 

LL'L" 

Therefore, we may obtain a stable result with the summations over a limited number of L when we consider the 
magnetic power spectrum is as red as ns ~ —2.9, because the summations of V and L" are limited by L. Here, we 
use the analytic formulas of the / symbols which are given by 

L" L 4 

£, ^ ^ 3 } oc (LL'L'T 1/6 , /C- 2 ( 2f +1) ocL"^ 1 , <j L 3 L 2 L : [. oc (L"L)-^ , (42) 



as described in detail in Appendix [B] 

Using the approximation and the summation rules described above, we can perform the computation of the CMB 
bispectrum containing full-angular dependence in a reasonable time. 



IV. RESULTS 

Now we show the result of the CMB temperature bispectrum induced from the vector anisotropic stress n^. 
In order to compute numerically, we insert Eq. (|38|) into the Boltzmann code for anisotropies in the microwave 
background (CAMB) [l2|, HH. We use the transfer function of magnetic-compensated modes calculated as Refs. 
[l4l l39i]. which is shown in Appendix [A] In the calculation of the Wigner-3j, 6j and 9j symbols, we use a common 
mathematical library called SLATEC |4(| and analytical expressions in Appendix [Bl 

In Fig. [2 we show the reduced bispectra of the temperature fluctuation induced by the PMFs defined as (25| 

Avvv) , T o \-i R (vw) ,,„v 

for 4 = 4 =4- Red solid and green dashed lines correspond to the bispectrum given by Eq. (1551) with the spectral 
index of the power spectrum of PMFs fixed as jib = —2.9 and —2.8, respectively. One can see that the peak of each 
bispectrum is located at £ ~ 1500 and the position is similar to that of the angular power spectrum Cj) induced 
from the vector mode as calculated in Appendix [C] At small scales, the vector mode contributes to the CMB power 
spectrum through the Dopplcr effect. Thus we can easily find that through this Dopplcr effect the vector mode can 
also enhance the CMB bispectrum. In our related paper [24j, we have also shown that the contribution from the 
vector mode dominates over that from the scalar mode at the scales around £ ~ 1500 in the CMB bispectrum induced 
from the PMFs, as in the CMB power spectrum. 

As for the amplitude of the CMB bispectrum of the vector mode induced from the PMFs, one can expect b^^m ~ 

( V)3 /2 

Cj p by using the amplitude of the CMB power spectrum of the vector mode induced from the PMFs. However, in 

Fig. [3 we find that the amplitude of b^J^ul ^ s smaller than the above expectation. This is because the configuration 
of multipoles, corresponding to the angles of wave number vectors, is limited to the conditions placed by the Wigner 
symbols. 

We can understand this by considering the scaling relation with respect to £. If the magnetic power spectrum given 
by Eq. (0| is close to the scale-invariant shape, the configuration that satisfies L ~ L" ~ I and L' ~ 1 contributes 
dominantly in the summations. Furthermore, the other multipoles are evaluated as 



Li ~ L[ ~ L'{ ~ I , L 2 ~ L'2 ~ L 3 ~ L' 3 ~ £ , L 2 ~ L 3 ~ 1 , 



(44) 
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FIG. 2: (color online). Absolute values of the normalized reduced bispectra of temperature fluctuation for a configuration 
l\ = li = £3. The lines correspond to the spectra generated from vector anisotropic stress for tib = —2.9 (red solid line) and 
—2.8 (green dashed line), and primordial non-Gaussianity with /Jj ^ = 5 (blue dotted line). The strength of PMFs is fixed to 
-BiMpc = 4.7nG and the other cosmological parameters are fixed to the mean values limited from WMAP-7yr data reported in 
Ref. 0. 



from the triangle conditions described in Appendix [Bl Then we can find b^^Y,} oc ^ 2 ™b+ 4 for £ < 1000, where we 
have also used the following relations 



k' 2 dkTj^(k)j t Mr -n))oci, { % **, £ l jeer 1 , IC-^ +1) ~ /C-^ +1) oc £ 

r 3/2 





n B +l 



(45) 



which, except for the first relation, are also coming from the triangle conditions of the Wigner 3-j symbols. Therefore, 
combining with the scaling relation of the CMB power spectrum mentioned in Appendix O we find that b^XY^J is 

suppressed by a factor £(™ B_1 )/ 2 from cj^ 3 ^ 2 . This is the reason why our constraint on the PMF from the vector 
bispectrum is not so much stronger than expected from the scalar counterpart. 

From this figure, we also find that the CMB bispectrum becomes steeper if ub becomes larger, which is similar to 
the case of the power spectrum. This will lead to another constraint on the strength of the PMFs. In particular, as 
shown in Refs. |18l - l22l . |24| . although the CMB bispectrum induced from the PMFs is dominated by the contribution 
from the scalar mode on large scales, such contribution becomes small on small scales. Therefore, it will be important 
to consider not only the contribution from the scalar mode induced from the PMFs on large scales but also that from 
the vector mode on small scales to obtain the constraint on the amplitude and the spectral index of the PMFs' power 
spectrum simultaneously. 

In Fig. [31 we show the reduced bispectrum bm^e^ with respect to £3 with setting i\ = i%. From this figure we 
can see that the normalized reduced bispectrum of the vector mode induced from the PMFs for ^1,^2^3 ^ 100 is 
nearly flat and given as 



4(4 + i)4(4 + l)|^;,:U|~2xio 



(WW I ~, o v m-19 



4.7nG J 



(46) 



It is seen that b^JJ^^ for tib — —3 dominates in i\ = £2 ^> ^3- This means that the shape of the CMB bispectrum 
generated from the vector anisotropic stress of the PMF is close to the so-called local-type configuration if the power 
spectrum of the PMF is nearly scale invariant. We can understand this by the analytical evaluation as follows. As 
mentioned above, in the summations of Eq. (|38p. the configuration that L ~ l\,L' ~ 1 and L" ~ £3 contributes 
dominantly. By using this and the approximations that 

L x ~ h , L[ ~ l 2 , L'l ~ h , L 2 ~L' 3 ~L, L' 2 ~L'l~L' , L' 2 ' ~ L 3 ~ L" , (47) 
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FIG. 3: (color online). Absolute values of the normalized reduced bispectra of temperature fluctuation given by Eq. (|38[) and 
generated by primordial non-Gaussianity given by Eq. (|48[) as a function of £3 with £\ and £2 fixed to some value as indicated. 
Each parameter is fixed to the same value defined in Fig. [21 



which again come from the triangle conditions from the Wigner symbols, the scaling relation of £3 at large scale is 
evaluated as &/j/X<2< x ^3 B+1 - From this estimation we can find that £i(£i + 1)^3(^3 + 1)&/j/X] 2 < x ^3 'S f° r 
tib = —2.9, and i® 2 for ns = —2.8, respectively, which match the behaviors of the bispectra in Fig. [3J 

In order to obtain a valid constraint on the magnitude of the PMF, we compare the bispectrum induced from the 
PMF with that from the local-type primordial non-Gaussianity in the curvature perturbations, which is typically 
estimated as [42[ 

i x {h + 1)£ 3 (£ 3 + l)b lll2l3 ~4x 10- 18 /NL Cal ■ (48) 

By comparing this with Eq. (l46)l . the relation between the magnitudes of the PMF with the nearly scale-invariant 
power spectrum and /^l^ is derived as 

(l^f) ~ 774 l/ ^ al|1/6 (fOT 11B ~ " 3) ' (49) 

By making use of the above relation, we can place the upper bound of strength of the PMF. If we assume | /nl al I < 
100 as considered in Ref. [lj|, we can translate this to the constraint on the PMF amplitude as BiMpc < 17nG, 
which is stronger by a factor of 2 than estimated in Ref. 19]. On the other hand, from the current observational 
lower bound from WMAP 7-yr data mentioned in Sec. Ill namely f^i, 3 * > — 10, we derive -Bim P c < UnG. If we use 
I /nl al I < 5 which is expected from Planck experiment [23(, we will meet a tight constraint as -BiMpc < 10nG. 



V. SUMMARY AND DISCUSSION 



In this paper, we present a calculation method of the bispectrum of CMB temperature fluctuation induced from 
the vector mode of the PMFs as described in Ref. 33], by taking into account the full angular dependence of 
the bispectrum of magnetic fields. We expand all the angular dependence with the spin spherical harmonics and 
convert them to the summations of the Wigner symbols. In the radial integrals and timelike integrals, we use only 
one approximation that the interval of the timelike integrals is confined to the moment of the recombination, which 
corresponds to neglecting the vector mode ISW effect. This approximation is valid because the radiation transfer 
function of the vector magnetic mode has a sharp peak around k ~ £ / (to — t*) which comes from the Dopplcr effect 
of the baryon vorticity induced from the magnetic field. We checked that the errors by the approximation are less 
than 10%. 

As the results, it is found that the CMB bispectrum from the magnetic vector mode dominates at small scales 
compared to that from the magnetic scalar mode which has been calculated in the literature. It is also found that 
the bispectrum has significant signals on the squeezed limit, namely the local- type configuration, if the magnetic field 
power spectrum is nearly scale invariant. This is understood by considering the asymptotic scaling relation of the 
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CMB bispectrum. We also investigate the dependence of the spectral index of the power spectrum of the PMFs on 
the CMB bispectrum and we find that the CMB bispectrum of the vector mode induced from the PMFs is more 
sensitive to the spectral index of the PMFs' power spectrum than that of the scalar mode. Hence, we conclude that 
it is important to consider not only the contribution from the scalar mode of the PMFs on large scales, but also that 
from the vector mode on small scales to obtain the constraint on the amplitude and the spectral index of the PMFs' 
power spectrum simultaneously. 

By translating the current bound on the local-type non-Gaussianity from the CMB bispectrum into the bound 
on the amplitude of the magnetic fields, we obtain a new limit: -BiMpc < HnG. This is a rough estimate and a 
tighter constraint is expected if one considers the full I contribution by using an appropriate estimator of the CMB 
bispectrum induced from the primordial magnetic fields. 

Because of the complicated discussions and mathematical manipulations, here we restrict our attention to the 
temperature bispectrum from the vector mode of the PMFs. However, one will be able to apply this methodology to 
the bispectra of CMB temperature and polarization from the scalar, vector and tensor modes. 
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Appendix A: CMB temperature fluctuations induced from vector anisotropic stresses 

In Refs. [H, H^, it is discussed that the temperature fluctuations are generated via Doppler and integrated Sachs- 
Wolfe effects on the CMB vector modes. Based on them, we derive the transfer function of the vector magnetic mode 
as follows. 

When we decompose the metric perturbations into vector components as 

Sgoc = Sg c0 = a 2 A c , (Al) 

Sg cd = a 2 {d c h { p +d d h^P) , (A2) 

we can construct two gauge- invariant variables, namely a vector perturbation of the extrinsic curvature and a vorticity, 
as 

V = A - h' , (A3) 
n = v A , (A4) 

where v is the spatial part of the four-velocity perturbation of a stationary fluid element and a dash denotes a partial 
derivative of the conformal time r. Here, choosing a gauge as h' = 0, we can express the Einstein equation 



V ' + 2 -V = - "^■ uv " 7 2 " v ' > , (A5) 

and the Euler equations for photons and baryons 

rJ 7 ' + T^(v 7 -v b ) = , (A6) 



«„' + ^fi b - ^(v 7 - v b ) = ^° n ^ V) . . (A7) 
a R a 4 {p b +p b ) 



Here = — iki,P ac Ti bc , p is the isotropic pressure, the indices 7, v and b denote the photon, neutrino and baryon, 
t c is the optical depth, and R = (pb + Pb)/(p-y +P-y)- In the tight-coupling limit as v 7 ~ v b , the photon vorticity is 
comparable to the baryon one: f2 7 ~ ri b = O- Then, the Euler equations (|A6I) and (IA7[) are combined into 

(l + R)Q' + R-n= fc ^ oII B V) , (A8) 
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and this solution is given by 

n(k,r) ~ /3(fc,r)n^ V) (k) , (A9) 



for k < ks 



P(k,r) = { (1+ 5 t^; ° +P - 0) , , / , (A10) 



fc(P"Y,0+P-y.(>) 



for k > ks 



where ks means the Silk damping scale. 

As mentioned above, the CMB temperature anisotropics of vector modes are produced through Doppler and inte- 
grated Sachs- Wolfe effect as 



AT(fi) 



h\H + J ° drV • n , (All) 



where tq is today and r* is the recombination epoch in conformal time, [ih.n = k • n, x = k(ro — r), and n is an unit 
vector along the line-of-sight direction. Because of compensation of the anisotropic stresses, a solution of the Einstein 
equation (|A5|) expresses the decaying signature as V oc a~ 2 after neutrino decoupling. Therefore, in an integrated 
Sachs- Wolfe effect term, the contribution around the recombination epoch is dominant. Furthermore, neglecting 
dipole contribution due to v today, we can form the coefficient of anisotropics as 



a im = J d 2 h^^Y; m (h) 



(2tt) 3 

In the transformation n — > (fJ,k,n, <pk,n), the functions are rewritten as 



< V) (k) • n -> -isj 1 —^ £ ng(k)e^ , (A12) 

V 1 A=±l 

y; m (n) D ^ (s&j) Y; m ,(n k<n ) , (ais) 

m' 

d 2 h -> dn k , n , (Ai4) 

where we use the relation: n^ y) = £ -ffl^e^ 3 and the Wi gner D matrix under the rotational transformation of 

A=±l 

an unit vector parallel to z axis into k corresponding to Eq. (A7) of Ref. [33| . Therefore, performing the integration 
over Q,k : n in the same manner as Ref. [HJ 3 , we can obtain the explicit form of ai m and express the radiation transfer 
function introduced in Eq. (|13j) as 



(£+i)j 



V2 



This is consistent with the results presented in Refs. [l2|, |3£ 



Appendix B: Analytic expressions of the Wigner symbols 

The Wigner-3j, 6j and 9j symbols exp ress Clebsch-Gordan coefficients between two other eigenstates coupled to two 
three, and four individual momenta [3#j36j . Their selection rules and several properties are reviewed in Ref. [28l[37|. 
Here, using their knowledge, we show the analytical formulas of the Wigner symbols which appear in the CMB 
bispectrum of Eq. (|3"8]). 



3 In Ref. [33l| . there are three typos: right-hand sides of Eqs. (B21), (B22) and (B23) must be multiplied by a factor — 1, respectively. 



The I symbols, which arc defined as iffff*' = y^ +i)(2'2+i)(2'3+ i| f ^ k ] . alv ( , xl)) ,, ssr( i as 



L hhh 



jO 1 -1 

/i h h 



(Eti l i)W(-h + h + h)W(h ~h + h)W(h +I2- l 3 y. 
(^±k). (Ai^+k)! (^^)! V /(EL 1 ^ + 1)! 



(for h+h + h = even) 



(for h+h + h = odd) 



' (k-l)(k + l) 
h - 1/2 

Z 2 + 1/2 



16tt 1 ' V (Z 2 - 1/2)(Z 2 + 3/2) 



(for ; x = h - 2, Z 3 = 2) 
(for/ 1 = ? 2 ,i 3 = 2) 



/± (_!)'» Mh + V ( io rh=h + 2,h = 2) 
8tt V ; V Z 2 + 3/2 V 1 2 ' 3 ' 



'-(-l) ls+ Vi 3 + l (for Ii =/ 3 - M2 = 1) 

<37T 



' — (-l) ,s+ Vis + l/2 (for « x = Z 3 ,Z 2 = 1) 



(-1)' 3+1 V^ (forli = l 3 + l,/ a = l) 



Three Wigner-9j symbols in Eq. (1380 are calculated as 



ti h h 
h h h 
1 1 2 



2(Z 3 ± 1) + 1 J h k 1 

; 3 ±2 j 3 ±i k 



h h 1 

ia±l ^3 h 



(for Z 6 = l 3 ± 2) 



' (2Z 3 - l)(2/ 3 + 2)(2Z 3 + 3) f h h 1 
30(2Z 3 )(2Z 3 + 1) \h h-l h 



/2(2Z 3 - l)(2/ 3 + 1)(2Z 3 + 3) \ h h 1 

^3 ^3 ^5 



15(2Z 3 )(2Z 3 + 2) 



h h 1 

J3-I h h 



h h 1 
^3 ^3 h 



+ \ 



'(2Z 3 - l)(2Z 3 )(2Z 3 + 3) f h h 1 



30(2; 3 + 1)(2Z 3 + 2) \h h + 1 h ) { h + 1 h h J (f ° r k ~ h) 



h h 1 
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where these Wigner-6j symbols are analytically given by 



(h h n = ( _ 1)tl+u+ , a+1 h+u+ i 6+ 2P2 i 1+h -i e+ iP2 {foxh = h _ lik = h+l) 

L '4 «5 «6 J y 2/4+3 "3 2/1+I "3 

Wx+U+ia+i +/4 + /e + 2)(/i +/ 4 ~/e + Zi + Z4 + k + l)(li - *4 + M 



(-1)' 
(-1) 



2/4+3^3 2/1+2-P3 

(for l 2 = l u l 5 = U + l) 



-1 -h+h+h + lP? /i-/ 4 +/6 + 1^2 



(for l 2 = h + l,h = h + 1) 



y 2U+3P3 2/1+3-P3 

/ 1+ /4+/ 6 + l [; 4 (; 4 + + h{h _ l){h + J) _ ; 6 (; 6 + X) _ ^ + ^ 



2(/i + Z 4 + + l)(/i + k - /e) 



" VHl+^+^6 + l)(/l-/4 + y2/4+2P3 2/ 1 +lP3 ^ h ~ h ^ ~ ^ 

2 ( 1 )' 1 +U+/ fi +i kfc + 1) + + 1)(Z 4 + 1) - Z 6 (Z 6 + 1) - Zi(*i + ljjj 

\/ 2/4+2-P3 2/1+2-P3 

(for J 2 = Z 1( Z 5 = '4) 
( _ 1) / 1 +/4+/ 6 + l ^ + !) + (/l + l){h + 2)(/4 + !) _ ^ + ^ _ ^ + 1)L J 



2(-h + h + h)(h -h + k + l) 



(for Z 2 = h + l,h = h) 



(Bll) 

(B12) 
(B13) 

(B14) 

(B15) 
(B16) 



y (Zi + z 4 + k + 2)(ii + z 4 - z 6 + 1)2/4+2^3 2/ 1+3 f:3 

Using these analytical formulas, one can reduce the time cost involved with calculating the bispectrum of Eq. (|38|) . 



Appendix C: CMB all-sky power spectrum of vector modes generated from PMFs 

In this section, we derive CMB power spectrum of vector modes sourced from PMFs in the same manner as presented 
previously and check the validity of our original approach. 

From Eq. (fl3|) , the CMB power spectrum of the intensity mode induced from n^ 1 ^ is formulated as 



a [V) a {V) * \ 



.71=1 



(2nf 



Ai,A 2 =±l 



Therefore, we should simplify the initial power spectrum of II^ 1 ^ as 



n 



ko rkn 

'^AV D_ (V \ I V^AV V_(V '\ I Al-Cl I A^-Cl 



(Cl) 



< 1 k^dk'P B (k[) J k' 2 2 dk' 2 P B {k' 2 ) / d 2 V x / cn4<5(ka - k' x - k 2 )£(k 2 - k' 2 - k^) 
x^ia4~ Al) (ki)fe c ei A2) (k 2 ) 



(C2) 
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For the first part in two permutations, we calculate <5-functions and the summations with respect to a, &, c, d: 

/•OO 

5(k 1 -k' 1 -k' 2 ) = 8 / A 2 dA (-l)^ 1± ^^lll% 3 JLAhA)jUk' 1 A)3LMA) 

JO T..T..T.. 



S(k 2 -k> 2 -K) = 8j a B dB YJ ^^lIvM^ B ^l'M B )H^ B ) 

M[M' 2 M's (C3) 

^^^^.^(-^(S -M^ -mJ > 
k la e d ^\k 2 )P ad {k\) = YJ E (y) ^ m S^\)x 2 Y lm A^)-^ ma {k\) a Y^ md {k\) , 

fc 2 4- Ai) (k 1 )n c (k^)= yj E (y) (-Ai^^-M^Cki)-^^^)^^^), 

ct'=±1 m c ,m b =±l,0 ^ ' 

perform the angular integrals of the spin spherical harmonics: 

[ d 2 k' Y* Yr » f Y* Y*. - Y^ (--[\ a+ma jO-v-S jO-a-S ( L' 3 1 L \ ( L 2 1 L \ 

f rl 2 V V* V V* V* — f— l^o-'+mc rO-a'-S' rO-a'-S' ( L 3 1 L' \ ( L' 2 1 L' 

jaV^WW^V^-^Ll ^ j l 3 il' ^ — M 3 m6 M' y y M 2 —m c M' 

/j2]T \r y T^* T^* \ , rOAl— Sk T" 

« > t\-\ 1 I\m b Ilm a -\i I i imi I l, 1 Mi — 2-j L 1 l 1 L k 1 \ 



L x l x L k \f 1 1 L, 
V Afi "2i Mfe M m m 6 M fe 

L k M k S k 



(C4) 



^ 2 l/ V V V V* — C _1 \™2 + A 2 tOA 2 -S p T 0X 2 -S p ( L[ l 2 L p \ I 1 1 Z/p 

LpA ly, Sp 



sum up the Wigner-3j symbols over the azimuthal quantum numbers: 



E (-D 



M 2 +m a 



MiM 2 M 3 
Mkm a mb 



1 1 L fc \ / Lx L 2 L a \ f L a 1 L' \ ( L 2 1 L\(L 1 l 1 L k 
m a m b M k ) \ Mi -M 2 -M 3 i I -M 3 m 6 M' / I M 2 -m„ M ) \ Ah m x M k 



I M / _ M mi M 1,3 ia L i 



L' L ^i 

J3 L 2 -^1 

1 1 Lk 



E 



M( M' 2 



+mc f 1 1 L p \fL[ L' 2 L' 3 \f L' 2 1 L' \( L> 3 1 L W f 2 L p 

-m c -m d M p i I M{ -M' 2 -M' z J I M' 2 -m c M' J I -M' z m d M J I -M[ m 2 M v 



C5) 



A/'+<? 2 +£o+L + l+L p f ^' ^ ^2 

M' -M m 2 





L 


/'2 




L's 


L\ 


1 


1 


L v 



and sum up the Wigner-3j symbols over M,M': 



MM' 



\rn-2 



E (-1) M+M ' f m' -M mi ) ( M — M m 2 ) = ^^i^A^.^ ■ ( C6 ) 



Following the same procedures in the other permutation and calculating the summation over L p as 

l + f-ll^ f L' L i 2 ) „ f V L i 2 \ 

E^i^^^V^ ^ ^ g - -^ 7 °& A2 L 'i L 'i L i ' (C7) 
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we can obtain the exact solution of Eq. (IC2[) as 



T (>2) 



2tt f 8(2n)^ 2 

3p 7 ,0 



*E E 

LL L1L2L3 



L, 



/(2^ + 1)<^ 

3 ^i+^-j 



-r 



jO 





L 










(? 


L 2 


£11 


!? 








1 






1 





^dAj Ll (fciA) / B A dB 3v ^k 2 B) 



' fc?d^P B (^)i £a (^>l)j L ;(fc^) / ° k' 2 2 dk' 2 P B (k' 2 ) jL ,(k' 2 B) jL3 (k' 2 A) 

Jo 

E(_t\L 2 +L' +L 3 + L' 3 rOS-S rOS-S T 0S'-S' tOS'-S' 
V J--* J L' 1L 1 L 2 1L 1 L 3 1L< 1 L' IL' 



S,S'=±1 



OA1-A1 7-OA1-A1 tOA 2 -A 2 



XMMl LieiLk l 11Lk l Vi 



\i 2 i 



(C8) 



In the 



Note that in this equation, the dependence on the azimuthal quantum number is included only in 5 mit 
similar discussion of the CMB bispectrum, this implies that the CMB vector-mode power spectrum generated from 
the magnetized anisotropic stresses is rotationally-invariant if the PMFs satisfy the statistical isotropy as Eq. 
Furthermore, using such evaluations as 



E (-D 



L 2 +L' 2 +L 3 + L' 3 jOS-S rOS-S tOS'-S' tOS'-S' 
1 L'1L 1 L 2 \L 1 L 3 1L' 1 L'1L> 



S,S'=±1 



<\iL€fAlib^\ib (for L> 3 + L 2 , L 3 + L> 2 = even) 
(otherwise) 

r01-l /-01-1 r01-l 



r 0A 2 -A 2 
\l 2 2 







= even) 



(otherwise) 



ErOAi-Ai rOAi-Ai tOA 
I L 1 £ 1 L k 1 llL k 1 L' 1 

,A a =±l 

2 poo 1 2 7t roo poo 

II 4 W ifrt^likn) / A 2 dAjL 1 (kiA) / B 2 dBj L{ (k 2 B) 

x / fci 2 ^P B (fci)j L2 (fciA)j L ,(fcii?) / k%dk 2 P B (k' 2 )3 L > 2 (k 2 B)j L3 (k' 2 A) 
Jo Jo 



(C9) 
(CIO) 



II 4 - 



( 27r ) 



xA|(r -r,) 4 (^ 



-(n B + l), 



(Cll) 



the CMB angle-averaged power spectrum is formulated as 

2 
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00 1.2 



11 



(2^)3 



T* \ 2 (ns+1). 



x E E E E A Uro r.)' ) /C-^ +1) (ro - r.JJC^fo, - n 



L2L2 
L3L3 

























1 




1 





19 



10000 




FIG. 4: (color online). The CMB power spectra of temperature fluctuation. The lines correspond to the spectra generated 
from vector anisotropic stress of PMFs as Eq. (|C12[) (red solid line) and primordial curvature perturbations (blue dotted line). 
The green dashed line express the asymptotic power of the red solid one. The PMF parameters are fixed to nj = —2.9 and 
-BiMpc = 4.7nG, and the other cosmological parameters are fixed to the mean values limited from WMAP-7yr data reported in 
Ref. [H|. 



This has nonzero value in the configurations: 

[L k , Lx) - (2, \£ ± 2|), (2,1), (1, £), L[ = \l ± , 
\£-L\ < L' < £ + L , 

(L 2 ,L' 3 ) = (\L-l\,\L±l\),(L,L),(L + l,\L±l\) , 

(L' 2 ,L 3 ) = (\L'-l\,\L'±l\),(L',L%(L' + l,\L'±l\), 1 1 

L\ + L 2 + L3 = even , L[ + L' 2 + L' 3 = even , 

\L! -L 2 \<L Z <L X +L 2 , \L[ - L' 2 \ < L' 3 < L[ + L' 2 . 

This shape is described in Fig. |4j From this figure, we confirm that the amplitude and the overall behavior of the 
red solid line are in broad agreement with the previous studies (e.g. [H, [H, EE El ) • For £ < 2000, using the scaling 
relations of the Wigner symbols at the dominant configuration L ~ £, L' ~ 1 as discussed in Sec. IIV1 we analytically 
find that oc £ nB+3 . This traces our numerical results as shown by the green dashed line. 
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